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requirements  for  the  Master  of  Science  degree  in  Naval 
Architecture  and  Marine  Engineering  and  the  Professional 
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ABSTRACT 

The  problem  of  applying  the  theory  of  systems  analysis 
to  the  marine  propulsion  plant  is  investigated  in  this  thesis. 
It  is  believed  by  the  author  that  the  use  of  this  theory  can 
be  of  considerable  assistance  to  the  naval  engineer  in  giving 
him  an  understanding  of  the  fundamental  behavior  of  the  system 
and  will  provide  a  scientific  basis  for  design  criteria  and 
operational  policies. 

There  are  five  basic  areas  of  investigation  covered  in 
the  thesis  which  all  contribute  to  the  organization  of  the 
overall  method  of  accomplishing  an  operational  analysis  of  a 
complex  mechanical  system.  The  following  is  a  concise  digest 
of  each  of  the  sections. 

1,  The  derivations  of  the  basic  mathematical  formulas 
for  reliability  determination  are  given.  The  concept  of  com- 
ponent failure  distribution  is  presented  and  the  exponential, 
normal,  and  Weibull  failure  models  are  discussed, 

2,  The  importance  of  defining  the  objectives  of  the 
analysis  clearly  and  completely  is  shown.  The  interpretation 
of  component  operating  conditions  and  the  resolving  of 
success  and  failure  standards  are  explained  and  an  outline 
for  the  formulation  of  the  necessary  input  data  required  to 
undertake  the  mathematical  part  of  the  analysis  is  offered, 

A  hypothetical  combined  diesel-gas  turbine  propulsion  plant 
is  used  to  illustrate  the  proper  procedures  for  performing  a 
reliability  analysis. 

3,  The  mathematical  procedures  for  determining 
system  reliability  by  the  correlation  of  the  individual 
system  components  behavior  are  explained.  In  addition  to 
presenting  methods  for  solving  serial,  parallel,  and  Bayesian 
systems,  formulas  are  derived  for  the  reliability  analysis  of 
periodically  operating  components  and  of  standby  units. 
Representative  values  are  postulated  for  the  diesel-gas 
turbine  plant  and  reliability  curves  are  plotted  which  demon- 
strate the  effect  of  component  wear-out,  reliability  as  a 
function  of  time,  and  the  individual  effect  of  specific  com- 
ponents on  the  overall  system  reliability. 
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4 .  Two  techniques  for  ascertaining  maintenance 
schedules  are  formulated.  The  first  technique  yields  the 
optimum  time  to  perform  preventive  maintenance  in  order  to 
maximize  system  availability  or  minimize  expected  repair 
costs.  The  second  procedure  deduces  by  dynamic  programming 
an  optimal  decision  policy  on  whether  or  not  to  repair  the 
system  dependent  upon  the  number  of  trips  made  since  the  last 
repair  and  the  number  of  trips  still  to  be  accomplished. 

5 .  A  brief  investigation  is  made  into  the  problem 
of  components  subject  to  the  interaction  of  their  failure 
distributions.  Four  possible  interaction  functions  of 
instant  failure  rate  are  proposed  and  analyzed  as  to  their 
effect  on  reliability  and  mean  time  to  failure.  An  example 
is  used  to  show  how  important  this  interaction  effect  can  be 
when  it  is  time  dependent, such  as  would  be  caused  by  gradual 
wear. 

In  conclusion  this  thesis  presents  a  procedure  for  ana- 
lyzing a  complex  mechanical  system  using  hypothetical 
examples  to  better  illustrate  the  utility  of  the  results  to 
the  naval  engineer. 


Thesis  Supervisor:   Ernst  G.  Frankel 

Title;  Assistant  Professor  of  Naval  Architecture 
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NOMENCLATURE 

Meaning 
Reliability 
Instant  failure  rate 
Probability  density  function  of 

failure 
P(t)>  Q(t)  Probability  distribution  of 

failure 
Tf ,  MTBF  Mean  time  before  failure 

t3   t  Time 

\  Exponential  instant  failure  rate 

C  Constant 

Specific  nomenclature  for  diesel-gas  turbine  reliability  analysis 
C.,  p(t)  for  diesel  engine 

Cp  p(t)  for  gas  turbine 

Co  p(t)  for  heat  exchange 

Cu  p(t)  for  reduction  gear 

Cc  p(t)  for  lube  oil  pump 

Cg  p(t)  for  standby  lube  oil  pump 

C7  p(t)  for  fresh  water  pump 

Co  p(t)  for  salt  water  pump 

N  Mean  time  of  non  use  of  the  gas 

turbine 
Cp  Probability  that  the  gas  turbine 

will  not  be  in  use  at  time  t 
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INTRODUCTION 

Operational  Analysis  and  Reliability  Engineering 

During  and  shortly  after  the  Second  World  War*  due  to 
the  growing  complexity  of  all  types  of  operations  in 
general,  a  branch  of  mathematics  termed  operations  analysis 
became  prominent.  Probability  theory  forms  the  basis  for 
this  field  of  science.  If  a  definition  is  desired*  it  may 
be  said  the  operations  analysis  is  a  scientific  study  of 
operations  to  better  understand  their  behavior  for  the 
twofold  purpose  of  predicting  their  future  results  due  to 
changes  in  the  system  and  controlling  the  operation  to  im- 
prove its  result.  This  thesis  is  concerned  with  the  use  of 
operations  analysis  methods  to  solve  problems  in  the  field 
of  reliability  engineering. 

In  any  system*  mechanical  or  otherwise*  one  of  the  most 
important  parameters  determining  its  overall  value  is  the 
reliability  that  the  system  can  be  expected  to  have  during 
its  operation.  Although  the  exact  definition  of  reliability 
is  usually  dependent  upon  who  is  the  def iner*  the  accepted ' 
definition  endorsed  by  the  American  Society  for  Quality 
Control  [  i  ]  states  that  reliability  is  the  probability  of  a 
device  operating  within  specified  limits  for  the  time  and 
operating  conditions  imposed  on  the  device'.  Operations 
analysis  techniques  are  used  in  reliability  engineering  for 
the  purpose  of  predicting  reliability  and  optimizing  the 
system  from  a  reliability  standpoint.  While  reliability 
analysis  is  a  relatively  newa  subject  which  has  generally 
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only  been  applied  to  the  electronics  and  missile  field,  it 
has  the  potential  of  being  extremely  useful  in  the  under- 
standing and  improving  of  all  types  of  mechanical  systems 
including  marine  propulsion  plants. 

Importance  of  Reliability  Analysis 

The  ability  to  be  able  to  predict  reliability  is  tanta- 
mount to  understanding  the  underlying  causes  of  unrelia- 
bility„  The  consequences  of  unreliability  are  numerous |  the 
most  obvious  ones  being  cost,  lost  time,  and  danger  to 
personnel o  The  cost  involved  not  only  is  that  of  the  device 
to  be  replaced  or  repaired,  but  also  the  cost  of  the  mainte- 
nance in  material,  manpower,  and  training  with  the  added 
expense  of  having  to  continuously  monitor  the  equipment  for 
possible  failures „  The  value  of  the  analysis  is  that  in 
choosing  some  specific  system  to  perform  a  task,  a  relia- 
bility analysis  of  the  system  will  furnish  data  on  which  to 
base  a  decision »  For  instance,  the  Armed  Services  mainly 
places  the  emphasis  for  selecting  a  component  or  system  on 
the  low  bid  price „  Normally,  reliability  is  only  given 
token  consideration  in  the  decision „  The  question  arises  as 
to  whether  a  reliability  investigation  of  the  device  should 
be  made  since  isn't  it  possible  that  it  would  be  advisable 
to  pay  more  initially  for  a  reliable  system  that  in  the  long 
run  will  give  better  availability  and  less  maintenance  costs 
than  a  cheap  unreliable  system?  Conversely,  the  possibility 
arises  as  to  whether  in  some  applications  too  much  relia- 
bility is  required  of  a  device .  Certainly  the  adage  that 


"a  chain  is  only  as  strong  as  its  weakest  link"  applies  in 
a  system  and*  for  example*  designing  a  valve  with  a  safety 
factor  of  5  when  the  piping  only  has  a  safety  factor  of  2 
is  a  waste  of  time  and  money.  It  is  therefore  evident  that 

reliability  apportionment  should  enter  into  any  preliminary 
system  design. 

Another  problem  inherent  to  this  day  and  age  is  that  of 
complexity o  The  major  problems  to  be  solved  have  become 
more  intricate  as  time  passes  and  it  is  only  natural  that 
there  is  an  increase  in  the  complexity  of  the  equipment 
designed  to  meet  these  problems .  The  more  complex  equipment 
is *  the  greater  number  of  components  are  required,  and  the 
more  likely  it  is  that  one  of  these  components  will  fail. 
Thus*  an  increase  in  complexity  leads  to  a  decrease  in 
reliability.  Of  what  value  is  complex  equipment  with  theo- 
retically high  performance  if  failures  of  its  components 
keep  the  piece  under  repair  most  of  the  time.  It  can 
therefore  be  seen  that  engineering  improvements  for  relia- 
bility will  in  many  cases  improve  ©perability*  maintaina- 
bility and  productivity  with  an  overall  net  gain  in  system 


Other  System  Analysis  Results 

Results  other  than  reliability  may  also  be  obtained  from 
an  operations  analysis  of  a  system.  The  question  of  main- 
tainability of  the  system  should  normally  be  investigated 
unless  for  some  special  reason  the  system  cannot  be  main- 
tained during  its  operating  cycle*  i.e.  missiles*  torpedoes* 
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ete.  A  study  of  maintainability  will  give  answers  to  such 
questions  asj 

1.  Is  preventive  maintenance  worthwhile  for  the 

particular  system? 

2,  If  worthwhile 3   at  what  time  intervals  should  it  be 

scheduled? 

3»  What  is  the  optimum  efficiency  that  can  be  expected 
from  a  system  preventively  maintained? 

In  conjunction  with  preventive  maintenance  scheduling* 
the  renewal  theory  of  operations  analysis  [  2.  ]  will  give  the 
expected  number  of  replacements  for  a  component  over  its 
lifetime^  t„  This  gives  a  good  basis  for  making  an  estimate 
of  the  number  of  spare  parts  to  be  stocked,, 

Some  other  results  that  are  obtainable  arej 

1,  The  effect  of  redundancy  on  reliability . 

2.  The  system  availability  that  can  be  expected. 

3o  Optimization  procedures  for  reasons  of  reliability 
or  coste 

40  Confidence  levels  that  can  be  applied  to  the 
results  of  the  analysis. 

5o  Information  for  the  determination  of  equipment 
quality  control  specifications. 

It  is  seen  that  by  the  methods  of  operations  analysis, 
a  deep  insight  into  the  behavior  of  a  system  can  be  obtained 
which  will  produce  not  only  a  better  understanding  of  the 
system  but  also  give  a  firm  basis  for  making  decisions  to 
improve  its  overall  value. 
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Present  State  of  the  Art 

At  the  present  time  the  state  of  the  art  of  operational 
analysis  of  mechanical  systems  is  very  sketchy.  Many  mathe- 
matical models  have  been  devised  of  hypothetical  systems 
which  have  been  analyzed  thoroughly  but  the  analysis  of  real 
systems  has  only  been  accomplished  in  the  electronics  and 
missile  fields  and  even  there  only  in  limited  situations , 
There  are  two  main  reasons  for  this  lack  of  progress.  First 
and  foremost  is  the  fact  that  the  value  of  systems  analysis 
and  especially  reliability  engineering  has  just  in  the  past 
few  years  come  into  prominence.  Now  that  the  mathematical 
tools  have  been  developed  and  the  value  realized,  the  majori- 
ty of  large  companies  are  beginning  to  organize  reliability 
analysis  departments.  The  second  factor  and  possibly  the 
hardest  to  solve  is  the  cost  and  time  involved  in  running 
the  testing  programs  necessary  for  the  analysis.   In  order 
to  obtain  any  results  from  an  analysis ,  failure  tests  must 
be  run  for  the  devices  under  investigation  and  reams  of  data 
must  be  evaluated.  Although  this  statistical  aspect  of 
operations  analysis  will  not  be  discussed  in  this  thesis s 
it  should  be  mentioned  that  there  are  sampling  procedures 
that  can  be  used  effectively  to  make  data  correlation  ef- 
ficient. Inherently.,  the  cost  of  initiating  a  reliability 
program  will  be  high  but  just  as  inherently  Is  the  fact  that 
in  the  long  run,  the  program  will  more  than  pay  for  itself 
in  its  results. 
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Thesis  Objective 

The  object  of  this  thesis  is  to  show  how  the  methods  of 
operations  analysis  can  be  employed  in  the  investigation  of 
a  marine  propulsion  plant  and  of  what  worth  such  an  investi- 
gation will  be.  The  overall  study  of  this  question  will  be 
divided  into  five  parts j 

1.  Reliability  analysis  basic  concepts. 

2.  Formulation  of  the  problem  of  determining  system 
behavior  by  correlation  of  individual  system  component 
behavior. 

3*  Mathematical  determination  of  system  reliability. 

4.  Effect  of  preventive  maintenance. 

5.  The  problem  of  interacting  components. 
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PROCEDURE 

I.  Reliability  Analysis  Basic  Concepts 
Reliability  and  Ms  an  Time  to  Failure 

Before  introducing  the  method  by  which  to  solve  the 
reliability  analysis  of  a  complex  propulsion  system,  a 
short  explanation  of  the  theory  behind  the  mathematical 
derivation  of  reliability  will  be  given.  This  is  deemed 
necessary  since  reliability  analysis  is  essentially  a  new 
field  and  although  it  has  as  a  foundation  the  mathematics 
of  probability  and  of  statistics,  there  are  some  intrinsic 
differences.  The  standard  reliability  terms  are  defined  in 
the  Nomenclature  section  of  this  thesis,  however  since  the 
basic  mathematical  formulation  is  short,  it  will  be  included 
in  the  text. 

To  determine  reliability  mathematically,  the  following 
is  presented; 

Let 

x  =  f (x, ,x0, „ . „,x)  be  a  vector  of  performance 
x     d.  n 

characteristics  of  the 
component 
A  =  region  of  the  sample  space  representing  the 

satisfactory  performance  of  the  component, 
w  =  f  (w-pWo* .  o  .,w  K  be  vectors  specifying  the 

>  upper  and  lower  limits  on 
z  =  f{z1,z2,  <>oo,zn)  J  satisfactory  performance 
t  =  time  of  the  interval  under  investigation 
f(x,t)  =  probability  density  function  of  the  performance 
characteristics 
R  =  reliabilitv  of  the  comoonent 
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therefore 

R(x  C  Aj  t)  =  Prob.(w  <_  x  <  z)  =  /  f(x,t)  dx 

A 
employing  the  proper  definitions  [zo] 

R(t)  =  /°°  f(x)  dx  =  1  -  /*  f(x)  dx  =  1  -  F(t) 

t.  0 

furthermore       _;t  p(x)  ^ 
R(t)  =  e  ° 

where 

p(t)  =  §{11=  -  ft  (ln  R(t))o 

-/  p(x)  dx 
f (t)  =  p(t)e  ° 

The  other  parameter  of  reliability  mathematics  that  is 

of  great  importance  is  the  components '  mean  time  before 

failure  (MTBF) „ 

If       Tf =  mean  time  before  failure 

E  =  expected  value 

then 

T     =  E(t)   =  /°°  t  f (t)    dt  =  J°°  t   dH(t)   =  /°°  R(t)    dt 
1  0  0  0 

In  the  event  that  the  component  has  been  operating  for  a 

period  previous  to  the  time  of  investigation; 

Defining  t  =  0  time  component  is  new5  or  system  starts 

operating 

t  =  t-,  time  of  investigation 

t  =    time  to  point  of  interest 

R(t/t, )  =  reliability  of  the  component  at  time  t  given 

the  component  has  operated  to  time  t^. 

thus 


R(t/tn)    =■   [t 


15' 


and 


R(t 


Tf   =   /       (t-t-^    flfT  f (t)    dt 


=  skr  f  R(t)  dt 

Component  Failure  Distributions 

From  the  preceding  paragraphs  It  can  "be  seen  that  relia- 
bility and  mean  time  to  failure  can  be  evaluated  from  the 
probability  density  life  distribution  of  the  component.   In 
the  following  paragraphs  some  of  the  more  Important  mathe- 
matical approximations  of  component  life  distributions  are 
discussed.  These  distributions  have  been  derived  from  ex- 
tensive analysis  of  operating  data  and  life  tests.  It  must 
be  remembered,  however*  that  they  are  just  statistical  models 
characterizing  a  physical  phenomenon  and  are,  at  least,  a 
good  approximation  of  the  reality.  For  a  better  understanding 
of  the  real  situation,  confidence  limits  are  associated  with 
each  of  these  approximations. 

Normally,  the  characteristic  parameter  used  to  describe 
life  distributions  is  the  instant  failure  rate,   p(t).  The 
time  dependence  of  p(t)  that  is  representative  of  a  normal 
component  is  shown  in  Fig.  I. 
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FIGURE  I 
Representative  p(t)  vs  t  plot 


.ft) 


where 


I  =  early  failure  region  which  results  from  manufacturing 
flaws s   in-transit  damage,  etc.  which  cause  normal 
failure  at  an  early  stage  in  the  operating  cycle. 
Usually  in  real  situations  there  is  a  "debugging"  or 
run-in  period  which  eliminates  this  problem. 
1^  =  chance  failure  region  caused  by  random  failures  which 
are  normally  considered  independent  of  time  and 
having  a  fixed  probability  of  oc curing  at  any  time  in 
this  region. 
Ill  =  wear-out  region  due  to  the  gradual  deterioration  of 
the  component  with  time. 

Exponential  Distribution 

This  is  the  most  widely  used  life  distribution  in 
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reliability  analyses „  The  reasons  for  this  are  it; 

1.  characterizes  components  that  have  been  optimized 
to  limit  failures , 

2.  characterizes  components  having  a  predominence  of 
human  errors, 

3.  characterizes  complex  components, 

4.  characterizes  components  consisting  of  parts  of 
mixed  ages, 

5.  approximates  what  is  usually  the  behavior  in  the 
chance  failure  region, 

6«  is  the  easiest  mathematically  to  work  with. 
The  equations  associated  with  the  exponential  distribution 
are  j 

p(t)  =  constant  ■  c 

f(t)  =  p(t)  e~P(t)t 

R(t)  =  e-P(t)t 

xf  pTtJ 

FIGURE  II 
Exponential  Distribution 


m 
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Normal  Distribution 

This  distribution  is  generally  employed  to  describe  the 
wear-out  phenomenon.  It  also  has  been  found  to  apply  to^ 
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1,  components  exhibiting  homogeneous  deterioration 
properties, 

2.  components  subjected  to  small  variations  in  environ- 
mental severity, 

3«  components  whose  failures  occur  at  times  well 
removed  from  t  =  0S      and  whose  MTBF  is  large 
compared  to  its  standard  deviation. 
The  equations  associated  with  the  normal  (Gaussian)  distri- 
bution are  i 

M  =  mean  lifetime 

a   =  standard  deviation  of  M 


't-M)2 

f  (t)   =  -i~    e         2a 

aV2? 


,2 


R(t)   =  -1—     Te         2<J         dt 
<W2lr     t 


Tf  =  M 

o  2 

if     erf (x)   =  ~  Jx  e~u     du 

v  0 

then 

R(t)  -  i  [1  -  erf  (£*) 


Normal  Distribution 
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FIGURE  III 


M     t 
Weibull  Distribution 

Although  this  distribution  has  rarely  been  used  in 

reliability  analyses  to  date,  statistical  investigation  of 

life  tests  of  mechanical  tests  shows  that  this  distribution 

seems  to  be  applicable  in  many  cases.  The  distribution  can 

be  described  by  a  two  parameter  method  or  a  three  parameter 

method o  It  includes  as  a  special  case  the  exponential 

distribution  and  can  approximate  a  normal  distribution 

(P  =  3.25)o 

For  the  three  parameter  casei 

where  a  =  scale  parameter 

P  =  slope  parameter 

y   =  location  parameter 


t  >  ys   a,    p  >  0 


P(t)  =  si* 


t  (t)  =  £&^ 
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(t-Y)P 

a 
R(t)  =  e 

T^  -  y  +  aP   r(p+l) 


Weibull  Distribution 


FIGURE  IV 


Other  Distributions 


t-v 


t-  v 


Some  other  distributions  occasionally  used  in  relia- 
bility studies  are j 

1 .  lognormal 

2.  gamma  or  beta 

3.  mixed  Weibull 

however ,  the  exponential 9   normal,,  and  Weibull  are  the  most 
important  of  the  distributions  and  are  the  ones  used  for 
examples  in  this  thesis. 
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II  Formulation  of  the  Problem 
Defining  the  Problem 

One  of  the  most  critical  parts  of  performing  an  analysis 
of  any  system  Is  the  complete  definition  of  the  problem  to  be 
studied.  This  definition  basically  must  answer  three 
questions | 

1.  What  are  the  objectives  of  the  analysis? 

2.  What  components  are  to  be  considered  comprising  the 
system? 

3.  What  will  be  considered  the  failure  events  of  the 
system  and  components? 

These  must  be  clearly  and  exactly  defined  if  clear  and 
exact  answers  are  to  be  expected.  In  some  cases  if  just 
figure  of  merit  or  trend  results  are  desired/  greater  liber- 
ties may  be  taken  in  the  definitions.  Once  these  three 
questions  are  answered,  the  problem  has  been  defined  and  the 
results  of  the  investigation  are  completely  dependent  upon 
these  definitions. 
Objectives 

When  defining  the  objective  of  the  analysis,  the  question 
is,  what  is  wanted  and  what  should  be  expected  from  the 
analysis.  In  order  to  better  illustrate  how  to  accomplish  an 
operations  analysis  of  a  system,  a  hypothetical  marine  pro- 
pulsion plant  will  be  postulated  and  investigated.  It  must 
be  realized  by  the  reader  that  this  example  will  be  necessarily 
simplified  since  the  object  is  not  numerical  results  but  . 
instead  the  method  of  the  solution. 
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The  objective  of  this  analysis  will  be  to  find  the 
reliability  of  a  combined  diesel  and  gas  turbine  marine 
propulsion  plant.  The  result  desired  is  how  does  the 
reliability  of  the  overall  plant  vary  with  time  in  oper- 
ation and  how  should  one  go  about  improving  this  relia- 
bility. Naturally  in  this  example  the  value  of  the 
numerical  results  are  useless  since  there  is  no  statisti- 
cal basis  for  the  numbers  used.  In  an  actual  analysis, 
statistical  life  tests  would  have  to  be  performed  on  all 
the  components  involved .  From  these  statistical  tests, 
confidence  limits  can  be  found  and  applied  to  the  final 
results.  Methods  of  finding  confidence  levels  are  dis- 
cussed in  Lloyd  and  Lipow  [ ?  ] . 
System 

The  next  step  in  the  problem  is  to  define  the  system. 
The  system  can  be  defined  in  as  basic  or  complex  terms  as 
the  investigator  deems  necessary  in  order  to  obtain  the 
results  desired.  In  this  example  the  system  will  be  com- 
prised of  the  following  components | 

1.  diesel  engine 

2.  lube  oil  pump 

3o  fresh  water  pump 

4.  salt  water  pump 

5.  heat  exchanger 
6C  gas  turbine 

7.  reduction  gear 

8.  standby  lube  oil  pump. 
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From  the  definition  of  reliability  it  is  seen  that  in 
defining  the  system*  the  operating  conditions  which  are 
within  the  prescribed  limits  must  also  be  enumerated. 
There  are  five  basic  operating  conditions  that  should  be 
taken  into  ac count j 

1 o  environmental 

2,  conditions  imposed  by  the  operator  or  user  of  the 
system 

3o  functional  dependence  of  the  components  within  the 
defined  system 

4.  interaction  of  component  failure  distributions 

5.  maintenance  policies. 

It  is  altogether  possible  that  more  than  one  combi- 
nation of  these  conditions  may  want  to  be  investigated.  In 
this  case  different  reliabilities  can  be  expected  for  the 
different  combinations.  For  the  purpose  of  this  investi- 
gation the  following  operating  conditions  are  defined^ 

Environmental s  The  only  extraordinary  environmental 
conditions  imposed  on  the  system  will  be  the  possi- 
bility of  shock  with  the  probability  of  P^  of 
occurring. 

Imposed  conditions ;  The  system  must  be  able  to 
operate  for  20$  of  its  total  operating  cycle  above 
cruising  power. 

Functional  dependence  i  The  functional  dependence  of 
system  operation  upon  the  components  isj 
1.  Diesel  must  perform  successfully  for  the  system  to 

operate . 
P   RArhir.-h-Inn  cpflr  must  nerform.  successfullv  for  the 
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system  to  operate. 
3o  Heat  exchanger  must  perform  successfully  for  the 
system  to  operate. 

4.  Salt  water  pump  must  nerform  successfully  for  the 
heat  exchanger  to  operate. 

5.  Fresh  water  pump  must  perform  successfully  for  the 
heat  exchanger  to  operate. 

6.  Gas  turbine  must  perform  successfully  when  more  than 
cruising  power  is  required. 

7«  Lube  oil  has  to  be  available  to  the  diesel  at  all 

times . 
8.  A  defined  shock  level  with  a  probability  of  P_  of 

occurring  will  cause  failure  of  the  reduction  gear 

and  heat  exchanger. 
Interaction  effects ;  No  interaction  among  the  component 
failure  distributions  will  be  assumed.  This  is  not 
necessarily  the  case  in  the  real  system  but  at  the 
present  state  of  the  art5  interaction  effects  cannot  be 
handled  except  by  using  "guess"  factors.  In  a  later 
part  of  the  thesis ,  these  interaction  effects  will  be 
investigated. 

Maintenance  policies :  For  the  propulsion  plant  of  this 
example.?  no  true  maintenance  actions  will  be  considered 
performed.  It  will  only  be  assumed  that  minor  adjust- 
ments for  power ,  temperature ,  and  fuel  regulation  will 
be  made.  Of  course  in  most  systems*  maintenance  actions 
will  improve  the  reliability  and  this  effect  will  be 
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covered  in  another  section  of  the  thesis. 
Concept  of  Success  and  Failure 

In  order  to  determine  the  failure  distribution  of  a 
component,  it  is  necessary  to  define  what  is  considered  a 
failure  and  what  a  success  in  the  operation  of  the  component. 
It  is  not  necessarily  true  that  the  only  failure  of  a  com- 
ponent is  when  it  stops  functioning  because  in  many 
situations,  the  system  will  also  fail  if  a  component  oper- 
ates outside  certain  tolerances.  For  instance,  if  a  lube 
oil  pump  doesn't  supply  the  oil  at  the  required  pressure 
even  though  the  pump  may  still  be  running,  the  system  will 
fail.  It  is  therefore  evident  that  the  investigator  must  be 
careful  in  choosing  the  ground  rules  for  success  or  failure 
of  any  test  and  relate  them  to  the  various  modes  of  operation. 

To  find  the  failure  distributions  relating  to  the  com- 
ponents, a  four  step  procedure  must  be  followed; 

1.  Define  success  and  failure  of  the  entire  system 
being  analyzed  for  a  particular  mode  of  operation. 

2.  Relate  the  definition  determined  from  1  to  each 
individual  components  operation  so  that  success -failure 
criteria  may  be  ascertained  for  each  component. 

3.  Establish  the  type  of  statistical  distribution  which 
describes  the  failure  phenomena. 

4.  Estimate  the  parameters  which  completely  define 
these  distributions. 

After  this  is  accomplished,  the  investigator  is  ready 
to  begin  the  mathematical  aspect  of  the  analysis  since  all 
input  data  are  now  available.  For  the  example  of  this  thesis, 
it  will  be  assumed  that  the  above  steps  have  been  carried  out 
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and  the  failure  distributions  of  the  components  in  question 
are  known.  The  case  where  the  out  of  tolerance  distribution 
is  independent  of  and  different  than  the  "catastrophic" 
failure  distribution  will  be  taken  into  account  in  one  case. 
Reliability  Schematic  Block  Diagram 

The  reliability  schematic  block  diagram  sometimes  called 
the  reliability  structure  model  is  a  picture  form  of  the 
functional  relationship  of  the  system  and  the  components 
under  investigation.   It  is  not  a  necessary  facet  of  an  oper- 
ational analysis,  however  in  complex  situations  the  model 
will  provide  an  insight  into  how  to  best  attack  the  problem 
from  the  mathematical  standpoint.  From  this  diagram  it  is 
also  easier  to  recognize  the  weaknesses  or  potential  weak- 
nesses of  the  system  from  a  reliability  sense.   The  relia- 
bility schematic  block  diagram  for  the  thesis  example  is 

shown  in  figure  V, . 

FIGURE  V 


*   DIESEL 


I  I 


B 


GAS  TURB. 


L.O.  PUMP 


STBYL.Q.P. 


F.W.  PUMP 


HEAT  EX. 
1 
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S.W.  PUMP 


RED. GEAR 


SHOCK 
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III  Mathematical  Procedure 

Three  Basic  Rules  of  Reliability  Mathematics 

Before  the  mathematical  analysis  of  the  system  is  begun, 
a  statement  of  the  basic  reliability  rules  of  probability 
calculus  will  be  discussed. 

The  first  rule  states  that  if  there  are  N  mutually 

independent  components,  each  having  a  reliability  of  Ri(t) 

and  that  for  satisfactory  operation  of  the  system  all  the 

components  must  function  properly;  then  the  overall  system 

reliability  R  (t)   is 
sys 

R  (t)  =  TT  Mt) 

sys    i=l   x 

This  is  known  as  the  Product  Rule  of  Reliability  and  a 
system  of  this  type  is  called  a  serial  system. 

The  second  rule  states  that  if  there  are  N  mutually 
independent  components  each  having  a  reliability  of  Ri(t)   and 
unreliability  of  Q±(t)   so  that 

R±(t)  +  Q±(t)  =  1 

and  that  for  satisfactory  operation  of  the  system,  only  one 

of  these  components  must  function  properly,  then  the  overall 

system  reliability  R  (t)  can  be  determined  since 

sys 

Q  (t)  =  TTQ^t)  =TT  (i-Mt)) 
sys    i=l  1      i=l 

R  (t)  =  1  -  Q  (t) 
sys        sys . 
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N 

r  (t)  =  i  -  TT  (l-Mt)) 

sys       i=l     x 

This  is  known  as  the  Product  Rule  of  Unreliabilities  and 
a  system  of  this  type  is  called  a  parallel  system.  This  is 
usually  described  by  the  term  redundancy^ relating  to  the  fact 
that  there  are  alternate  components  to  help  the  system 
operate  successfully  in  case  of  failure  of  one  or  more  of 
the  other  components.  The  normal  method  for  increasing 
system  reliability  is  to  add  redundancy  to  the  original  system. 

The  third  and  most  important  rule  is  based  on  Bayes ' 
probability  theorem.  It  states  that  the  probability  of 
system  failure  is  the  probability  of  system  failure  given 
that  the  component  is  bad  times  the  probability  the  component 
is  bad  plus  the  probability  of  system  failure  given  that  the 
component  is  good  times  the  probability  the  component  is  good. 

Qg(t)  =  Qg(t)|R1(t)  x  Rt(t)  +  Qg(t)|Q±(t)  x  Q±(t) 

This  is  called  the  Conditional  Rule  of  Reliability  and 
can  be  applied  to  any  system  including  serial  and  parallel. 
Mathematical  Analysis  of  the  Complex  System 

With  these  rules  the  analysis  of  the  complex  system  can 
be  begun.  They  will  not  account  for  every  situation  en- 
countered in  the  investigation  but  they  will  be  the  foundation 
for  all  the  computational  work  done. 

For  the  purpose  of  simplicity,  the  components  will  be 
referred  to  by  the  number  in  the  upper  left  hand  side  of  the 
block  in  Figure  V  ,   i.e.  1  =  diesel,  2  =  gas  turbine,  etc. 
It  is  seen  that  blocks  A,  Bs    G,   D  form  a  serial  system  such 
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that  if  any  one  of  these  fail  the  system  fails.  If  the 
reliability  of  each  of  these  four  blocks  could  be  found  then 

R  (t)  =  RA(t).RB(t).Rc(t).RD(t)  I 

This  will  be  the  approach  to  determ  ining  the  overall  system 
reliability. 
Reduction  Gear 

The  first  component  to  be  examined  will  be  the  reduction 
gear.  It  is  known  that  the  reliability  of  the  gear  is  Rj,  (t) 
under  normal  operating  conditions.  However,  it  is  also  known 
that  some  predefined  shock  amplitude  will  cause  the  gear  to 
fail. 

Define  P   as  the  probability  that  a  shock  of  the 
defined  amplitude  will  occur 

%   -  1  -   Ps 
Reliability  of  D  =  Reliability  of  4  given  no  shock 

occurs  times  the  probability  of  no 
shock  plus  the  reliability  of  4 
given  the  shock  occurs  times  the 
probability  that  shock  occurs. 
RD(t)  =  (R4(t)|PB  =  0)  Qs  +  (R4(t)|P,s)Ps 

(Vt)|l?s  =  0)  =  R4(t)   R4(t)|Ps  =  0 
therefore 

Rp(t)  =  R4(t)  Qs  II 

Heat  Exchanger 

In  finding  the  total  reliability  of  block  C,  we  again 
use  the  Bayes  Theory  approach  of  Rule  3. 
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Rc(t)(  Cs  =  0)=  R3(t)(if  7  is  good>R?(t)+R3(t)(if  7  is  bad)^(t) 

=  R3(t)(if  7  is  good)-R7(t) 

as  R3(t)(if  7  is  bad)  =  0 

which  leaves  R3(t)(if  7  is  good)  to  be  found 
R3(t)(if  7  is  good)  =  R3(t)(if  8  is  good)-Rg(t)+ 

R3(t)(if  8  is  bad)-Qg(t) 

R3(t)(if  7  is  good)  =  R3(t)(if  8  is  good>Rg(t) 

as  R3(t)(if  8  is  bad)  =  0 


but 


R3(t)(if  7  and  8  are  good)  =  Rc(t)(  Eg  =  o) 


,\  Rc(t)(Ps  =  0)=  R3(t>R7(t)-R8(t) 

and  it  is  seen  that  the  dependence  of  the  heat  exchanger  on 
the  fresh  water  and  salt  water  pumps  is  the  same  as  a 
subsystem  composed  of  these  three  components  in  a  serial 
arrangement.  However,  the  shock  hasn't  been  taken  into 
account  therefore 

Rc(t)  =  R3(t>R7(t)-Rg(t)-Qs  III 

Gas  Turbine 

In  the  analysis  of  the  block  containing  the  gas  turbine 
the  first  unusual  case  of  mathematical  reliability  prediction 
appears  where  the  three  basic  rules  are  not  sufficient 
enough  to  give  an  answer.  In  this  case,  probability  theory 
must  be  used  to  formulate  a  new  type  of  reliability  equation. 
This  is  due  to  the  fact  that  the  gas  turbine  is  only  used 
periodically  and  only  needs  to  operate  properly  when  called 
into  use* 
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The  reliability  of  the  gas  turbine  can  therefore  be 
found  by  the  following  procedure ; 
Define 

f2(t)  -  failure  density  function  of  the  gas  turbine 
P2(t)  =  probability  that  the  gas  turbine  is  not  in 

use  at  time  t_ 
U2(t)  =  distribution  function  of  periods  of  non-use  of 

the  gas  turbine 
Qo(t)  =  unreliability  of  the  gas  turbine. 

The  situation  can  be  described  by  the  statement^  a 
failure  of  the  gas  turbine  at  time  t  £  t  causes  system 
failure  by  t_    under  one  of  two  conditions: 

1.  The  gas  turbine  is  in  use  at  time  t  when  it  fails, 

2.  The  gasturbine  is  not  in  use  at  t  but  is  called 
into  use  after  t  and  at  or  before  t. 

Condition  ^L  is  characterized  by 

o 
Condition  2  is 

Qg(t)  =  J*  F2(T)U2(t-T)f2(T)dT 
and 

Qgtt)  =  Q2(t)  +  Qg(t)  IV 

Since  the  periods  of  non  use  are  considered  to  be 
randomly  spaced. 
Defining 

N2  =  mean  period  of  gas  turbine  non-use 
the  non-use  density  function  is 
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«  (t)        1  p"  ^ 
U2^)  =  !g  e 


and 


•    9* 


-  t/n, 


.   t     -  T/N0 

u2(t)  -/**-'•     dT  = 

o  1N2  o 

U0(t)  =  1  -  e"  t/N2 


2 
Diesel  Engine 

In  deriving  the  reliability  equation  for  block  A, 
encompassing  the  diesel  engine 3   another  situation  arises  when 
the  three  basic  reliability  equations  are  not  adequate.  This 
situation  is  caused  by  the  lube  oil-standby  lube  oil  pump 
combination  and  this  combination  is  referred  to  as  a  standby 
system.  To  find  the  reliability  of  such  a  system 
Define 

ff-(t)   =  failure  density  function  of  the  lube  oil  pump 

fg(t)  =  failure  density  function  of  the  emergency  lube 
oil  pump 
Since  the  standby  lube  oil  pump  only  operates  at  some  time 
t, 3     when  the  lube  oil  pump  fails*  the  true  failure  density 
functions  for  the  system  are 

fc-(t1)     lube  oil  pump 

fg(t-t.,)  =  standby  lube  oil  pump 

As  the  only  time  the  system  fails  is  when  both  pumps  fail, 

the  system  unreliability  is  found  by  Rule  2. 

Q  (t)   =  Qf-(t)Q6(t)  =  J*  f  (t^dt^  f6(T-t1)dx 
standby    D  o  J  o 

and  the  failure  distribution  function  of  the  emergency  lube 
oil  pump 
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^(t-t^  =  /*  f6(T-t1)dT 
i 

R(t)  =  1  -  Q(t) 
consequently 

R   (t)   =  1  -  J*  Fgft-tJMtJdt...  V 

standby      o   D    15  11 

Using  the  same  type  analysis  as  for  the  heat  exchanger 

RA(t)  =  ^(tjR  (t)  VI 

standby 

Overall  System  Reliability 

To  find  the  overall  system  reliability,  equation  I  is 

applied 


R  (t)  =  RA(t)*RB(t).Rc(t).RD(t)  1(a) 


where 


R.(t)  =  R,(t)*R  (t) 
R  x  standby 

RB(t)  =  [1  -  Q2(t)] 

Rc(t)  =  R3(t)-Rr(t)-R8(t).QB 

RD(t)  =R4(t).Qs 

Equation  I (a)  is  applicable  for  the  overall  system 
reliability  for  any  type  failure  distributions  &>r   the  com- 
ponents i.e.  exponential,  normal,  Weibull,  etc.  In  most 
instances,  however,  the  solution  of  such  complex  equations 
necessitates  the  use  of  a  digital  computer.  For  the  purpose 
of  easy  understanding  and  readily  obtainable  results,  the 
failure  distributions  of  the  components  will  be  assumed  expo- 
nential. This  is  the  standard  method  in  the  majority  of 
reliability  analysis  and  for  most  complex  systems  not  too 
unrealistic,  at  least  for  a  first  approximation.  The 
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assumption  of  the  exponential  distribution  requires  two  basic 
postulates; 

1.  The  components  are  not  operated  long  enough  so  that 
they  begin  to  wear  out  due  to  aging.  Prom  a  reliability 
standpoint  this  is  a  logical  way  to  operate  the  components 
as  their  reliability  greatly  decreases  as  they  enter  the 
wear  out  regime. 

2.  Any  failure  of  a  component  is  completely  random, 
i.e.  it  is  not  time  dependent. 

Using  the  exponential  distribution,  the  following 
equations  are  arrived  at; 

RA(t)  =  e   X  [e  5  +^-V  (e  5   -  e  6  )] 

6  5^ 


See  Appendix  [B  ]  for  proof 


t 


c  t  c      ~  W"~        c 

RB(t)  =  e"  2  [14-  P2(t)(1   2   )]  +  e   2  [P2(t)(r— 2— )] 

Nj  "  C2  1£  "  C2 

See  Appendix  [ B  ]  for  proof . 

Rc(t)  =  Qg[e   37  8] 

-Cj,t 

RD(t)  =  Qg[e  4  ] 

thus 
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t 

(   -Cpt    Pp(t)Cp       "  Np   Po(t)C9 
times  he  2  (l+  jS £)]  +  e    2  [^ £] 

I  N^  "  C2  N^  "  C2 

Quantitative  Analysis 

In  order  to  more  graphically  show  the  value  of  the  relia- 
bility analysis «,  representative  values  for  the  component 
instant  failure  rates   (po(t)  ~  c^)  will  be  chosen.  The 
overall  system  reliability  as  a  function  of  time  can  then  be 
evaluated. 


The  values  chosen  are; 

Cl  = 

.0002 

hr"1 

c2  = 

.0003 

hr"1 

C3  = 

.00006 

hr"1 

c4  = 

.00004 

hr"1 

C5  = 

.0004 

hr"1 

C6  = 

.0005 

hr"1 

C7  = 

.0005 

hr"1 

c8  = 

.0005 

hr 

Qs  = 

.98 

P(t)  = 

.80 

N2  = 

200  hr 

Also,  assume  the  wear-out  of  only  one  component  is  taken 
into  account.  This  component's  failure  due  to  wear-out  shall 
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be  described  by  a  normal  distribution  with  a  mean  life  time 

of  1200  hrs  and  a  varience  of  400  hrs. 

Figure  VI  shows  R  (t)  vs  t  for  the  cases  with  and  with- 

sys 
out  one  component  wear-out  consideration. 

For  the  purpose  of  understanding  the  individual  com- 
ponents relationship  to  the  system  as  a  whole,  each  of  four 
component  failure  distributions  have  been  varied  to  find  its 
total  effect  on  the  system.  These  four  components  were  chosen 
due  to  the  different  nature  of  their  operation; 

Component  Time  Operation 

Diesel  continuous,  dependent  upon  L.  0.  system 

Gas  turbine  periodic 

L.  0.  pump  continuous  with  standby  system 

F.  W.  pump  continuous,  nondependent 
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IV  Preventive  Maintenance : 

In  most  of  the  reliability  work  done  in  the  past,  little 
interest  has  been  given  to  the  effect  of  preventive  mainte- 
nance .  The  reason  for  this  is  twofold  due  to  the  fact  that 
most  of  the  reliability  study  has  been  applied  to  the  missile 
and  electronic  fields £ 

1,  In  the  missile  field  the  operation  is  usually  "one 
shot"  i.e.  the  component  or  system  having  once  begun  its 
operation    cannot  be  maintained,  and  after  terminating  this 
operation  it  is  unusable. 

2.  In  the  electronic  field,  the  failure  distributions 
of  the  components  are  considered  exponential  thus  making 
preventive  maintenance  useless.  This  is  so  because  the 
instant  failure  rate  is  independent  of  time  and  therefore 
replacing  the  component  at  some  specific  time  when  it  is 
still  operating  successfully  has  no  advantageous  effect  from 
a  reliability  standpoint. 

These  two  conditions  do  not  necessarily  hold  however^ 
for  a  marine  propulsion  plant.  Such  plants  definitely  are 
not  "one-mission"  devices  and  it  is  doubtful  if  the  exponen- 
tial distribution  will  be  applicable  to  the  mechanical  com- 
ponents found  in  such  systems.  It  is  proven  by  R.  F.  Drenick 
["=?  ]  that  if  the  instant  failure  rate  is  constant  or  decreases 
with  time,  then  preventive  maintenance  is  not  beneficial.  In 
mechanical  components,  it.  has  been  determined  that  the  instant 
failure  rate  normally  increases  with  time.  In  this  case  pre- 
ventive maintenance  can  be  worthwhile  and  should  be  considered 
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in  any  operational  analysis  of  a  marine  propulsion  system. 
A  maintenance  investigation  should  provide  two  basic 
parameters^  how  valuable  to  system  performance  is  preventive 
maintenance,  and  at  what  time  intervals  should  it  be  per- 
formed. For  the  purpose  of  analysis  two  different  replacement 
policies  will  be  explored.  The  mathematical  analysis  is  from 
the  work  of  R.  Barlow  and  L.  Hunter  [ 5  ] . 


After  every  t   hours  of  continuing  operation  without 
failure,  preventive  maintenance  is  performed.  If  there  is  a 
failure  before  t  ,   the  maintenance  is  performed  at  this 
time,   t, ,  and  the  next  maintenance  is  scheduled  at  time 

1    o 

Defining 

y\m    =  expected  fractional  amount  of  time  the  system 

is  operating  as  t  — >  °° 

the  two  parameters  that  determine  the  condition  of  the  problem 

are  r\m    and  t  .  The  criterion  of  optimality  will  be  the 

maximization  of  tj^  . 

T  =  expected  time  needed  to  perform  repair  to  the 
e 

system  after  a  failure 

T  =  expected  time  needed  to  perform  a  scheduled 
s 

preventive  maintenance  action. 
By  making  the  assumption  that  after  each  repair  action,  the 
system  is  as  good  as  new,  the  following  equations  are  appli- 
cable for  optimal  results » 
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p(t  )  A  R(t)dt  -  F(t  )  = 


o  e  s 


'I"  -  l+(Te-Tg;P(t0) 


e    s 


if  T  <  T  ,   no  preventive  maintenance  should  be  made  and 

6      S 

T 

^loo  "  T-+T 
i   e 

Policy  II : 

In  this  case,  preventive  maintenance  is  performed  on  the 
system  after  it  has  been  operating  for  a  period  of  t°  hours 
regardless  of  the  number  of  intervening  failures.  If  there 
is  a  failure  before  t°  hours  have  elapsed  only  minimal 
repair  is  made  which  in  effect  puts  the  system  back  in  oper- 
ating condition  but  does  not  change  the  basic  system  failure 
rate. 

Using  the  same  definitions  for  r\       and  T   and  defining! 

00       s 

T  =  expected  time  needed  to  perform  minimal  repair 
to  the  system 
the  following  equations  are  determined 

/   t  p»(t)dt  =  ^ 
o  m 

"S"  =  1+Tm  p(t°) 

The  question  naturally  arises  as  to  which  policy  should 
be  used  for  a  particular  situation.  This  question  can  be 
resolved  by  equating  ti1oo  =  r\2oo     for  a  given  Tg  which  is 
known  and  the  same  for  both  policies.  By  substitution  of 
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values  for  T  .  T   can  be  determined.  By  graphing 

v?      Ill 

T^  vs  T™  for  a  given  T,  which  policy  to  use  becomes 

6        III  S 

evident.  Figure  XI  shows  such  a  graph  for  a  Weibull 

distribution. 

A  Maintenance  Policy  by  Dynamic  Programming 

Assume  that  a  preventive  maintenance  policy  is  to  be 

determined  for  a  system  having  an  instant  failure  rate  that 

is  a  linear  function  of  its  operating  time.  The  system  is 

to  be  operated  in  discrete  time  intervals  of  length  t  . 

If  the  system  is  repaired  between  operating  periods  the  cost 

associated  with  the  repair  is  C  .  If  the  system  fails 

during  an  operating  interval  it  cost  CL  for  the  repair  and 

the  system  cannot  be  operated  again  till  the  next  interval. 

Define 

C  =  cost  of  scheduled  maintenance 
r 

Cf  =  cost  of  "in  service"  repair 

P  =  Probability  of  failure  in  interval  m  +  1 

m  ____„«_ 

given  it  has  not  failed  in  m  intervals 
since  the  last  maintenance 
f  (m)  =  least  expected  cost  of  making  n  more 
intervals  given  m  intervals  have  been 
made. 
There  are  two  policies  that  can  be  followed,  either 
repair  the  system  at  m  or  do  not  repair  it.  By  the 
methods  of  dynamic  programming  covered  by  R.  Bellman  [ 6  ] 
the  following  dynamic  programming  equation  is  developed^ 


4i. 


f  (m)  -  min 
n 


~R:  Cr+Po[Cf+fn-l(0)]  +  (l-Po)fn-l(l) '   repair 
NR:   P  [C-+f  ,  (0)]  +  (l-P.Jf  , (m+l);  non-repair 


m 

CL  >  C   for  logical  results  otherwise  preventive  mainte- 
nance is  not  worthwhile. 
To  show  mathematically  the  effect  of  the  failure  repair 
cost,   Cf,  and  the  scheduled  maintenance  cost,   C  ,  on  the 
maintenance  policy  to  be  followed,  the  following  problem  will 
be  postulated. 

A  system  has  a  failure  distribution  curve  that  can  be 
described  by  the  Weibull  approximation  with  the  following 
parameters  j 

y   =  location  parameter  =  0 

{3  =  shape  parameter  =  2 

a  m   scale  parameter  =  57  x  10  hrs 

a  '  P  =  characteristic  lifetime  =  r\   =  755  hr 
M(t)  =  N(t)  =  interval  time  =  100  hr 
therefore 

p(t) 


R(t)  =  e  57(10^) 

The  question  to  be  answered  is  how  does  the  maintenance 
policy  vary  as  a  function  of  m9   n,  Cf,  and  C"r„  Figures 
XII  through  XVI  give  the  results  of  this  investigation. 
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V  Interaction  of  Component  Failure  Distributions 

In  the  previous  sections  of  the  thesis,  the  overall 
reliability  of  a  hypothetical  complex  marine  propulsion 
plant  has  been  derived  by  the  proper  mathematical  combi- 
nation of  individual  component  reliabilities.  In  actual 
situations  however,  it  has  been  found  that  although  this 
method  is  mathematically  correct,  it  does  not  yield  the 
actual  reliability  observed  in  practice.  Intuitively,  it 
might  at  first  glance  appear  that  this  poor  correlation  is 
because  the  model  is  not  a  good  functional  representation 
of  the  real  system.  Although  this  is  a  possible  reason, 
further  analysis  may  show  that  mathematically  all  work  is 
accurate.  Where,  then,  is  the  mistake  in  the  analysis? 
One  of  the  major  possible  explanations  and  one  that  although 
recognized  is  not  included  in  present  reliability  studies, 
is  that  when  components  are  integrated  into  the  system, 
there  is  an  interaction  among  them  which  can  change  their 
failure  distribution.  Because  of  this,  the  failure  distri- 
bution observed  in  the  single  component  life  test  is  not 
applicable  and  must  be  modified.  The  problem  of  finding  what 
type  of  function  this  modification  factor  is,  will  in  actual 
practice  be  very  difficult.  However,  if  reliability  investi- 
gations are  to  be  accurate,  component  failure  rate  inter- 
action has  to  be  taken  into  account  as  it  may  have  a 
pronounced  effect  on  system  reliability,  preventive  mainte- 
nance policies,  and  system  availability  calculations.  It  is 
therefore  imperative  that  laboratory  life  tests  of  components 
suspected  of  interaction  and  data  assimilation  of  real  system 
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operation  be  analyzed  so  that  interaction  functions  can  be 
determined. 

For  the  purpose  of  illustration,  two  possible  situations 
where  interaction  may  play  a  substantial  role  in  an  investi- 
gation are  given. 

a.  The  failure  distribution  of  a  diesel  cylinder  liner 
will  be  affected  by  the  wear  of  the   piston  rings , 
the  wear  of  the  connecting  rod  bearings,  and  the 
fuel  quality. 

b.  The  failure  distribution  of  a  reduction  gear  will 

be  dependent  upon  the  wear  of  the  shaft  bearings  and 
coupling  operation. 
Since  no  data  is  available  on  interaction,  for  the 
purpose  of  investigation  only  hypothetical  interaction 
functions  can  be  postulated.  In  the  following  sections  such 
assumed  functions  are  examined. 
Effect  of  Interaction  on  Component  Reliability 

The  time  dependent  reliability  of  a  component  is  de- 
scribed by 

R(t)  =  exp[-  J*  p(t)  dT] 
o 

which  is  a  function  of  a  single  variable,   p(t)  =  the  instant 

failure  rate.  An  equation  illustrating  the  main  factors  that 

influence  component  failure  rate  is  characterized  by 

p(t)  =  f(a,  p,  <|),  t) 

a  =  derating  factor 

p  =  environmental  stress  factor 

<|)  =  interaction  factor 

t  =  operating  time 
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where 


a  =  accounts  for  the  fact  that  the  component  can  operate 

at  different  outputs 
P  =  is  a  function  of  the  operational  environment  the 

component  experiences,  i.e.  temperature,  vibration, 

and  acceleration  stress  levels 
4)  =  accounts  for  the  effect  of  interaction  among  the 

component  failure  distributions 
t  =  time  the  component  is  in  operation 

It  is  normally  possible  to  find  the  instant  failure  rate 
as  a  function  of  derating  factor,  stress  factor,  and  time 
from  experimental  life  testing  of  the  component.  An  example 
of  such  a  function  can  be  shown  in  a  diagram  for  a  component 
that  follows  the  exponential  failure  distribution. 
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Because  reliability  engineering  is  relatively  new,  little 
experimental  data  is  available  to  find  such  functional  re- 
lationships and  what  is  available  only  applies  to  the  expo- 
nential failure  distribution.  It  is  valid  to  assume  that  the 
same  type  of  relationships  or  equations  can  be  found  for  other 
types  of  components  and  failure  dsitributions .  The  problem 
still  presents  itself  •  how  can  the  interaction  factor  be 
handled?  Since  no  work  has  been  done  on  this  matter,  only 
conjectural  relationships  can  be  assumed  for  the  purposes  of 
investigation.  Such  relationships  will  therefore  be  assumed 
and  their  properties  analyzed. 
Define 

X(t)  =  component  instant  failure  rate  without  taking 
into  account  this  interaction  factor,  4). 
therefore 

p(t)  =  f  (  *(t),  (>) 

It  is  this  function  that  will  be  hypothesized.  In  all 
eases  x(t)  will  be  considered  a  constant  for  ease  of  mathe- 
matical computation.  The  effect  of  the  interaction  factor 
which  is  determined  in  this  exponential  case  will  apply  in 
effect  to  other  failure  distributions. 

Investigation  1:  <|)  =  additive  and  constant 

p(t)  =  A  +  2  C. 
n 

R(t)  =  exp[-(A  +  2  C.)t] 

n 

Tf  =   l/(*  +  2  C±) 
n 
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Investigation  2:  <|)  =  additive  and  time  dependent 

p(t)  =  \  +   2  C.t 
n 

R(t)  =  exp[-At]   exp[-  \  2  C±t2] 


n 


Tf  .  expt^L-]  [  7^-  (1  -  erf  [^_   -    I 

See  Appendix [Bj for  proof.  n 

Investigation  3 :  $   =  multiplicative  and  time  dependent 

p(t)  =  (ct)  A 

R(t)  =  exp[-  |  A  t2] 
1  2R4 

Investigation  k%     §  -   additive  and  a  function  of  a  constant 

and  a  time  dependent  term 
p(t)  =  *  +  C^  +  C2t 

from  correlation  with  investigation  2 

R(t)  =  exp[-  (VfC^t]  exp[-  \  C2t2] 


Tf  =  expt^#]j^(l-erf[^])} 


It  can  be  seen  that  these  values  for  R(t)  and  Tf 
can  differ  greatly  from  the  values  obtained  for  the  standard 
exponential  failure  distribution,  i.e. 

R(t)  =  exp[-At] 

Figure  XVII  shows  the  relative  importance  of  2  C.  for 

n 
investigation  2  where  \  -   .01. 
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Preventive  Maintenance  Policy  by 
Dynamic  Programming  for  C„  =  5 
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Preventive  Maintenance  Policy  by- 
Dynamic   Programming  for     Cf  =  10 

C     = 

2 

14 

13 

12 

n     =  ip 

11 

r         •* 

10 

■Id   V> 

9 

m 

8 

7 

6 

Repair 

5 

. 

1 

1 

4 

Non-Repair 

3 

w 

2 

• 

1 

r 

3 

i 

Cr  =  1//2                                         &&  4*2 

2          3         4          5 

6         7 

8          9         10 

n 

FIGURE  XIV 

-66- 

Preventive  Maintenance  Policy 

by 
15 

uynaimc 

JTA-U^,. 

l  dimiiJ-iig,    xvi         v~    — 

1  c 

1  c    = 

|    r 

15 

3 

14 
13 

12 

11 

10 

Q 

m 

8 

^p   =    e.    1/  [ 

• 

7 

6 

Vj 

r 

5 

4 

3 

Repair 

i 

l 

p 

r 

f 

<_ 

Non-Repair 

l 

•JU. 

0 

c„  =  1 
r 

j^ze  <?-6? 

] 

L 

2         3         4 

56789] 

L0 

n 

57' 


FIGURE  XV 

Preventive  Maintenance  roncy  oy 
Dynamic  Programming  for  Cf  =  20 
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Preventive  Maintenance  Policy  by 
Dynamic  Programming  for  ratio 
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DISCUSSION  OF  RESULTS 

Reliability  Analysis 

In  Figure  VI  the  overall  system  reliability  as  a  function 
of  time  is  shown.  The  plot  graphically  displays  what  degree 
of  reliability  can  be  expected  from  the  system  assuming; 

a.  the  values  for  the  instant  failure  rates  specified 
previously. 

b.  the  components  are  good  as  new  at  t  =  0. 

c.  no  preventive  maintenance  is  performed  on  the  system 
during  the  operating  cycle. 

It  can  be  determined  that  this  overall  reliability  curve 
for  the  system  consisting  of  different  components  with  unlike 
operating  conditions  can  be  approximated  by; 

R  (t)  =  e"Ct 
sys 

where        C  =  .00198 

Other  results  that  can  be  obtained  from  Figure  VI  are; 

1.  If  the  instant  failure  rates  of  all  the  components 

were  decreased  by  50$,  the  overall  system  reliability  will  be 

increased  greatly  as  t^  increases ,  i.e. 

t  =  200;  Re„0  is  increased  by  20.2$ 
sys 

t  =  500;  R  ro  is  increased  by  61.2% 
sys 

t  =  800;  R0,ro  is  increased  by  116.8$ 
sys 

2.  If  one  component  has  wear-out  characteristics  (fresh 
water  pump)  such  that  it  should  be  included  in  the  reliability 
analysis,  the  overall  system  reliability  is  adversely  affected. 
This  effect  becomes  especially  prominent  as  the  system 
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operating  time  approaches  the  mean  life  time  of  the  component 

and  is  also  dependent  upon  the  varience  of  the  mean  life  time. 

Figure  VI  shows  that  a*t  t  -  800,  R  __,  is  decreased  by  15.7$ 

sys 

due  to  wear-out.  The  equation  that  is  applicable  when  investi- 
gating a  system  or  component  which  exhibits  random  and  wear- 
out  failures  during  the  operating  cycle  isj 

R±(t)  =  R±  (t)-R±  (t) 
w     r 

defining  Rj(t)  =  overall  reliability  of  component  i 

R.  (t)  =  reliability  of  component  i  exhibiting 
w 

only  wear-out  failures 

R.  (t)  =  reliability  of  component  i  exhibiting 
r 

only  random  failures 

Since  the  wear-out  reliability  factor  is  multiplicative, 
the  system  reliability  will  decrease  greatly  with  time  as  the 
number  of  components  experiencing  wear-out  increases. 

3.  Figures  VII  through  X  illustrate  the  relative 
importance  of  each  component  to  overall  system  operation 
from  a  reliability  viewpoint.  These  figures  show  percentage- 
wise how  the  decrease  in  a  specific  components  instant  failure 
rate  will  affect  the  system  by  using  the  formulaj 

Rsys(wlth  decreased  p(t))-R(pi(t)=Ci) 
Rsys  =  Rsys^i(t,=Ci) 

It  is  observed  that  a  decrease  in  the  fresh  water  pumps 
failure  rate  by  50$  increases  the  overall  system  reliability 
by  22$  for  an  operating  cycle  of  800  hours  while  the  same 
decrease  for  the  diesel  engines  failure  rate  only  increases 
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the  reliability  by  7.8J6.  Therefore  if  the  system  reliability 
needs  to  be  improved,  strictly  from  a  reliability  standpoint 
of  the  four  components  investigated  the  fresh  water  pump  per- 
formance should  be  improved  first  and  then  the  lube  oil  pump, 
gas  turbine,  and  diesel  respectively.  An  investigation  like 
this  immediately  indicates  where  in  the  system  .improvements 
should  be  made  to  increase  overall  reliability.  Conversely, 
it  also  reveals  what  increase  can  be  expected  in  system  relia- 
bility for  an  improvement  in  a  specific  components  instant 
failure  rate.  Consequently,  graphs  of  the  type  shown  in 
Figures  VI  through  X,  which  can  be  determined  for  any  system 
if  the  instant  failure  rates  of  its  components  are  known, 
provide  an  intimate  insight  into  system  behavior  and  provide 
valuable  information  for  decisions  regarding  engineering 
improvements. 
Preventive  Maintenance 

Figure  XI  shows  the  relative  merit  of  using  either  of 
the  two  policies  discussed  on  pages  3q  -<qo       of  this 
thesis.  This  graph  is  representative  of  the  preventive 
maintenance  situation  for  a  system  having  a  failure  rate 
that  increases  with  time.  As  an  example  of  how  to  use  the 

figure,  assuming  T  =  7  hours  and  the  expected  value  of 

s 

T  =  4  hours  then  policy  I  is  applicable  unless  T  <_  2 
hours.  Similarly,  for  any  value  of  T   there  is  a  corre- 
sponding  value  of  T   which  determines  the  boundary  for 
choosing  the  optimum  policy  to  follow.  It  will  generally 
be  found  that  policy  II  is  most  applicable  to  complex  systems 
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which  have  many  subsystems  liable  to  fall  and  where  keeping 
records  of  times  to  failure  for  each  subsystem  becomes 
cumbersome.  It  is  interesting  to  note  that  if  costs  of 
failure  are  substituted  for  the  times  of  repair,  then  these 
same  equations  minimize  the  expected  cost  of  operating  the 
system. 

The  method  of  determining  preventive  maintenance 
scheduling  by  dynamic  programming  is  more  suitable  to  a 
marine  propulsion  plant  than  the  previously  discussed  method. 
The  reason  for  this  is  that  with  operating  schedules  to  main- 
tain it  is  probable  that  optimum  preventive  maintenance  times 
found  will  not  coincide  with  the  ships  schedule.  This  method 
also  provides  a  continuous  picture  of  what  policy  to  take  re- 
garding preventive  maintenance  over  all  combinations  of  time 
intervals  (trips)  completed  and  time  intervals  still  to  be 
completed.  Figures  XII  through  XVI  illustrate  the  results  of 
such  an  analysis  made  for  a  system  having  the  Weibull  failure 
distribution  described  in  the  Procedure  section.  From  these 
graphs  the  expected  cost  of  maintenance  can  be  minimized  for 
any  operational  situation  of  this  system  if  the  cost  of 
failure  and  of  maintenance  are  known.  These  figures  are  only 
applicable  for  m,  n  >  0  as  if  m  =  0,  the  system  has  just 
been  repaired  and  it  does  no  good  to  repair  it  again,  and  if 
n  =  0,  there  are  no  more  trips  to  make  so  the  expected  cost 

=  fi. 

Although  the  graphs  are  self-explanatory,  one  interesting 

observation  can  be  made  which  is  not  intuitively  obvious  nor 
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would  normally  be  expected  If  the  investigation  had  not  been 
made.  Take  as  an  example,  Figure  XIII  for  C~  =  10  and  the 
curve  representing  C  =  2.  It  is  seen  that  if  8  trips  have 
been  made  (m  =  8)  and  2  more  are  planned  (n  =  2),  the  policy 
should  be  to  repair  the  system.  However,  with  the  same  past 
history  and  if  3  more  trips  are  necessary  (n  =  3),  the  system 
should  not  be  repaired.  The  explanation  of  this  is; 

1.  at  large  m  the  probability  of  failure  in  n  inter- 
vals is  relatively  high  and  for  small  n,  repair  is  the  best 
policy  since  in  effect  it  starts  the  system  off  at  m  =  0. 

2.  as  n  increases,  this  effect  of  starting  at  m  =  0 
becomes  less  important  and  the  cost  of  preventive  maintenance 
is  the  dominant  factor.  The  reason  is  that  at  large  n,  even 
if  maintenance  is  performed,  there  is  still  a  significant 
probability  of  failure  in  the  following  n-1  intervals  which 
causes  a  greater  expected  cost  than  if  no  repair  is  made  and 

a  high  probability  of  system  failure  is  accepted  in  the  next 
n  intervals . 

Analysis  of  Figures  XII  through  XV"  provides  the  results 
for  the  determination  of  Figure  XVI.  Figure  XVI  shows  that 
policy  determination  is  solely  dependent  upon  failure  distri- 
bution, length  of  interval,  and  the  ratio  of  Cf  to  C  . 
This  fact  lessens  greatly  the  number  of  calculations  that 
have  to  be  made  in  a  system  analysis  since  assuming  the 
system  failure  distribution  and  interval  time  are  known,  only 
one  graph  has  to  be  derived  instead  of  one  for  each  possible 

value  of  C.p.  As  the  values  of  C«  and  C   are  bound  to 
i  i       r 
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vary  with  time  due  to  the  normal  trend  of  increasing  repair 
and  material  costs,  the  advantage  of  the  single  graph  repre- 
senting the  ratio  is  evident. 
Interacting  Components 

Although  little  if  any  investigation  has  been  made  into 
the  problem  of  the  interaction  component  failure  distri- 
butions, Figure  XVTI  illustrates  that  in  certain  cases  this 
effect  can  be  very  important.  For  instance,  it  is  seen  that 

if  p(t)  =  .01  +  4(l0""5)t  the  mean  time  to  failure  is 

Tf 
decreased  20$.  On  the  other  hand  for  t  =  ^—   ,  this  inter- 
action effect  only  decreases  reliability  by  less  than  1$. 
Consequently  although  the  effects  of  interaction  may  not 
affect  the  basic  calculations  for  reliability  significantly 
and  therefore  may  be  difficult  to  observe  and  determine,  they 
definitely  can  become  significant  in  the  overall  system 
analysis  through  Tf  which  is  used  in  many  operations  analy- 
sis calculations  not  discussed  in  this  thesis,  i.e.  renewal 
theory.  It  therefore  becomes  imperative  that  investigations 
should  be  made  determining  any  interaction  present  in  a 
system  if  worthwhile  results  are  to  be  obtained  from  an 
analysis  of  the  system. 
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CONCLUSIONS 

In  addition  to  the  specific  conclusions  stated  in  the 
Discussion  of  Results,  the  following  general  conclusions  can 
be  drawn. 

1.  The  operational  analysis  of  a  marine  propulsion 
plant  will  provide  results  not  otherwise  available  which 
will  be  beneficial  in  making  decisions  regarding; 

a.  preliminary  design 

b.  operational  policies 

c.  maintenance  actions 

d.  system  improvement  and  redesign 

2.  In  such  an  analysis  the  definition  and  formulation 
of  the  objectives  of  the  investigation  determine  all  suc- 
ceeding stages  of  the  analysis  and  the  utility  of  the  results, 

3.  Statistical  data  are  a  necessary  prerequisite  to  any 
operational  analysis.  The  definition  of  the  problem  governs 
the  areas  of  investigation  for  the  statistical  tests. 

4.  Reliability  as  a  function  of  time  can  be  character- 
ized by  the  single  parameter  p(t),   instant  failure  rate. 
The  effect  of  an  increase  or  decrease  in  this  parameter  is 
amplified  as  time  increases. 

5.  The  wear-out  failure  distributions  of  the  individu- 
al system  components  should  be  included  in  the  analysis  of 

a  mechanical  system  as  they  can  greatly  affect  the  overall 
system  reliability.  This  consequence  becomes  especially 
prominent  as  _t  approaches  the  MTBF  of  a  component. 

6.  From  an  operational  analysis  of  a  system,  the 
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relationship  of  each  component's  performance  to  overall 
system  operation  can  be  determined.  Such  information  will 
be  of  considerable  significance  in  the  determination  of 
policies  regarding  system  design  and  improvement. 

7 .  Preventive  maintenance  is  only  advantageous  for 
devices  having  an  instant  failure  rate  that  increases  with 
time. 

8.  For  this  type  device  there  is  an  optimal  time  for 
the  performance  of  preventive  maintenance  which  will  maxi- 
mize availability.  An  optimal  time  may  also  be  found  which 
will  minimize  expected  cost  of  operation. 

9.  Dynamic  programming  can  be  utilized  in  determining 
an  optimum  maintenance  policy  for  systems  operating  in 
discrete  time  intervals.  This  policy  will  give  a  complete 
picture  at  any  time  in  the  systems  operating  cycle  of  what 
maintenance  action  to  follow  in  order  to  minimize  the 
expected  cost  of  repair. 

10.  The  effect  of  interacting  component  failure  distri- 
butions normally  neglected  in  operational  analyses  and 
difficult  to  determine  quantitatively  can  considerably 
affect  the  results  of  the  investigation  and  should  therefore 
be  included. 
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RECOMMENDATIONS 

1.  The  concept  of  the  operational  analysis  of  marine 
propulsion  plants  should  be  accepted  by  naval  engineers  as 
an  important  tool  in  the  evaluation  of  existing  and  proposed 
marine  systems. 

2.  A  program  of  component  testing  and  data  assimilation 
should  be  initiated  immediately  for  the  purpose  of  determi- 
ning representative  failure  distributions  of  mechanical 
devices . 

3.  Accepting  the  fact  that  corroborating  statistical 
data  are  minimal,  systems  analyses  should  be  carried  out  on 
existing  systems  for  the  object  of  developing  and  refining 
techniques  and  gaining  facility  in  the  procedures  of  such  an 
investigation. 

4.  Further  examination  should  be  made  into  the  question 
of  how  out  of  tolerence  i.e.  wear-out  failures  are  to  be 
defined  and  handled. 

5.  A  detailed  investigation  should  be  made  into  the 
problem  of  failure  distribution  dependence.  This  should  be 
explored  both  from  a  theoretically  mathematical  and  an 
experimental  viewpoint. 

6.  The  applicability  of  the  method  of  dynamic  program- 
ming to  system  analyses  questions  should  be  investigated. 
Possible  applications  other  than  for  maintenance  policies 
might  be; 

a.  optimal  reliability  apportionment  among  the 
components , 
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b.  optimal  arrangement  of  components  under  con- 
straints for  maximum  system  reliability. 

c.  optimal  choice  of  components  to  minimize  the 
deleterious  effect  of  failure  distribution 
interaction. 
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APPENDIX  A 
Definitions 
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Failure :  a  detected  cessation  of  ability  to  perform  a  speci- 
fied function  or  functions  within  previously  established  limits. 
Independent  Failure :  those  component  failures  which  occur  or 
can  occur  without  being  related  to  the  malfunctioning  of  associ- 
ated components. 

Redundancy ;  the  existence  of  more  than  one  means  for  ac- 
complishing a  given  task  where  all  means  must  fail  before 
there  is  an  over-all  failure  to  the  system. 
Availability:  the  fraction  the  total  desired  operating  time 
that  a  system  is  actually  operable. 

Mean  Time  Before  Failure  (MTBF) :  the  mean  or  average  time 
between  successive  failures  of  a  component 

Cumulative  Probability  of  Failure  (Q(t),  F(t)):  probability 
that  the  life  of  a  component  is  less  than  t. 
Probability  Density  Function  of  Failure  (f (t) ) :  probability 
that  the  life  of  a  component  will  be  between  t  and  t  +  dt. 
Instant  Failure  Rate  (p(t)):  probability  that  a  component 
will  fail  during  the  time  interval  t  to  t  +  dt  conditional 
upon  it  surviving  up  to  time  t. 

Reliability  R(t):  probability  that  a  component  will  operate 
within  specified  limits  for  the  time,  t,  and  operating  con- 
ditions imposed. 
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APPENDIX  B 
Details  of  Procedure 

I  Periodically  used  component 

II  Standby  system 

III  Interaction 

IV  Computer  program 
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I.   Determination  of  R(t)   for  a  component  used  periodically 
assuming  an  exponential  failure  distribution. 

F(t)  =  J*  [l-P(T)]f(T)dT  +  J*  P(T)[0(t-T)]f(T)dT 

o  o 

P(t)  =  P 

(t-T) 

G(t-T)  =  1  -  e    N 
f(T)  -  C  e"CT 

Mt)  =  J*  [1-P]C  e"CTdT  -  [l-P][l-e-Ct] 
o 

=  1  -  P  -  e"Ct  +  P  e"Ct 

Q^t)  =  J*  P[l-e"¥     ]C  e"CTdT 
o 

t 

=  P  -  P  e_Ct[l  +  y-S— ]  +  ?(T°—)e~   ¥ 

N  "  C      N  "  C 

t 


F(t)  =  1  -  e-Ct[l  +  Hy^-J]   +  HT^)e"   N 
R(t)  -  1  -  P(t) 


IT  "  C      I  "  C 


R(t)  =  e-Ct[l  +  Hy^—  )]  -  Htt2—  )e"  * 

N  "  C       N  "  C 

II.  Determination  of  R(t)   for  a  standby  system  assuming 

exponential  failure  distributions. 

p(t)  =  C.,   for  main  component 
p(t)  =  Cp  for  standby  component 

For  a  standby  system  Q(t)  =  ^   F2(t-t1)f1(t1)dt1 
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*  +.   -cP(t-tn) 

F2(t-t1)  =  r   fft-t^dt  -  /*  C2e  *    -1  dt 

-C,  (t-t1 ) 

=  1  -  e  -1    -1 

-C  (t-t  )     -C  t 
Q(t)  =  /*  Fgft-^J^ftJdtj  =  /*  [1-e  2    1-]C1e  X  1dt1 
o  o 

o  o 

-C,t     C.    -CUt   (C~-a.)t-  t 
C2  Cl  o 

=  1  -  e  -1  -  TT-i-te  1     -  e  2  ] 
C2""°l 

R(t)  =  1  -  Q(t) 

R(t)  =  e  1   [1  +  E-4-]  -  e  2 

III.  Determination  of  R(t)  and  Tf  for  p(t)  =  A  +  Ojt 
p(t)  =  "K  +   C±t 

R(t)  =  expf-/1^*  +  C±T)dT] 
o 

=  expt-/*  XdT  -  /t  C.t  d-r] 
o        o 

=  exp[-*t  -  C±   J*  t  dx] 

o 

ci    p 

=  exp[-*t  -  £=■  t  ] 

1  r  i-2 
=  e    e 


»oo 


Tf  =  /  R(t)dt 
o 
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=  /°°  exp[-(M;  +  |  C±t2)]dt 


o 
From  C.R.C.  Standard  Mathematical  Tables  [7  ] 

o 
it  can  be  seen  that 

|  c.t2  +  u  -  I  c±[t2  +  f£  t(+  £-  2  -  £-  2)] 

i      >  2   >2 


therefore 

2  2 

R(t)  =  exp[^r-]  exp[-  |  C±(t  +  $-)  ] 

therefore 

2  2 

Tf  =  /°°  R(t)dt  =  exp[^-]    f°°  exp[-  |  C±(t  +  $-)    ]dt 

let  x  =  t  +  ^-  t  =  0  x  =  ■£ 

dx  =  dt  t=»  x  =  °° 

J00  exp[-  |  C±(t  +  ^-)   ]dt  =  /°°  exp[-  |  C±(x2)]dx 
o  i  \ 

which  equals  (dropping  the  subscript) 

r  exp[§  x2]dx  -  f*  exp[|  x2]dx 
o  o 

from  C.R.C.  tables 

J00  exp[§x2]dx  =^ 


now  let         t  =  v?y 
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=  ^  x  x  =  0  t  =  0 


dt  =  Jl  dx  x  -  \  t  -  -^ 


therefore 


j°  «p[§  x2]<te  -  £-><£>  4    I™    exp[-t2]dt 

o  vTf-  O 


O 


and  ^ 


S--^    exp[-t2]dt  =  erffJy 

sTf  o  ^ 

thus 

f exp[§  x2)dx  =  4g  [1  -  erf  (-^)] 


o 
and 


Tf  =  exp[^]     4j  (1  -  erf  {-£* ) 
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XV*  Main  Body  of  Computer  Program  for  Solution  of  Equation 
(FORTRAN) 

DO  114     N  =  1,  Nl 
JJ  =  0 

DO  114     M  =  1,  Ml 
IF(N-l)   100,  100,  101 

100  Al  =  0. 
GO  TO  102 

101  J  =  N  -  1 

Al  =  C0STl(j,l) 

102  A2  =  CASREP  +  Al 
A3  =  RATE(M)  *  A2 
IF(N-l)  103,  103,  104 

103  A4  =  0. 
GO  TO  105 

104  K  =  M  +  1 

A4  =  C0STl(j,K) 

105  A5  =  l.-RATE(M) 
A6  =  A5  *  A4 
C0STl(N,M)  =  A3  +  A6 
IF(N-l)  106,  106,  107 

106  Bl  -  0. 
GO  TO  108 

107  Bl  =  C0ST2(J,1) 

108  B2  =  SKDREP  +  RATE(l)  *  (CASREP  +  Bl) 
IF(N-l)  109,  109,  110 

109  B3  =  0. 
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GO  TO  111 

110  B3  =  C0ST2(j,2) 

111  B4  =  l.-RATE(l) 
B5  =  B4  *  B3 
C0ST2(N,M)  =  B5  +  B2 
IP(JJ)  112,  112,  113 

112  JJ  =  1 
PRINT  20 
PRINT  21,  N 

113  L  =  M  -  1 

114  PRINT  22,  L,  C0ST1(N,M),  L,  C0ST2(N,M) 
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